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ABSTRACT 

Ice  drift  in  the  marginal  ice  zone  (MIZ) 
is  a  very  important  feature  of  air-ice-ocean 
interaction  at  high  latitude.  Thermally  gener¬ 
ated  surface  winds,  blowing  from  ice  to  water 
(ice  breeze)  with  some  deflection  due  to  the 
earth  rotation,  force  the  ice  drift  and  ocean 
currents  near  the  MIZ.  By  changing  the  surface 
temperature  gradient,  the  ice  motion  and  the 
ocean  currents  feed  back  on  the  surface  winds. 
A  coupled  air-ice-ocean  theoretical  model  for 
the  MIZ  is  employed  to  discuss  the  ice  drift 
pattern  with  such  a  feedback  mechanism.  The 
steady-state  solutions  show  that  an  off-ice 
and  divergent  wind  field  not  only  produces  a 
dilation  of  the  MIZ  (as  people  generally 
think),  but  also  generates  a  compaction  of  MIZ 
for  some  circumstances.  An  ice  divergence/con¬ 
vergence  criterion  is  found.  The  time-depen- 
dent  solutions  show  that  the  ice  motion  exhib¬ 
its  two  bifurcations.  First,  it  bifurcates 
into  decaying  and  growing  modes.  Second,  the 
growing  mode  bifurcates  into  non-oscillatory 
and  oscillatory  states.  Finally,  the  model 
predicts  the  ice  edge  upwelling. 

NOMENCLATURE 

ice  concentration 
integral  constant 
integral  constant 

air  drag  coefficient  on  ice,  0.0036 
air  drag  coefficient  on  water, 

0. 0012 

water  drag  coefficient  on  ice, 

0. 00086  m/s 

mean  surface  temperature  difference 
across  the  MIZ  , 

gravitational  acceleration,  9. 81m/s* 
reduced  gravitational  acceleration, 

0. 005m/s- 

mean  ice  thickness 

equilibrium  depth  of  the  pycnocline 
thickness  of  the  upper  ocean 
wave  number 
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L  length  scale  twice  the  MIZ  width, 

200  km 

R  rate  for  Rayleigh  friction  and  Newton¬ 

ian  cooling,  2.08xl0“5/s 
Ri  atmospheric  Richardson  number 

U  scale  of  ice  breeze  wind 

ua,va  ice  breeze  wind 

u^,Vi  ice  drift  velocity 

u^v^^  ocean  cuurent 

u  angle  between  surface  wind  and  surface 

temperature  gradient 

P  angle  between  ice  velocity  and  surface 

wind 

6  thickness  of  atmospheric  Ekman  layer 

C  nondimensional  ice  edge  displacement 

0a  surface  temperature 

0_n  characteristic  surface  temperature, 

270°K 

X^-i  eigenvalue  for  vertical  structure 

eigenvalue  for  time  evolution 
v  air  vertical  eddy  viscosity,  5  m*/s 

pa  air  density,  1.29  kg/m3 

ice  density,  920  kg/m3 
p„  water  density,  1000  kg/m' 

i l  angular  velocity  of  the  Earth  rotation, 

0.  7292x10“  vs. 

INTRODUCTION 

The  importance  of  ice  breeze  in  the  MIZ 
is  discussed  as  follows.  Schmidt  (1947)  used 
a  simple  linear  model  to  depict  seabreeze/ 
landbreeze  phenomena  and  concluded  that  for  a 
maximum  land-water  temperature  gradient  of 
4.  3"K/100km,  a  maximum  landbreeze  intensity  of 
about  2  m/s  would  be  reached.  From  observa¬ 
tions  in  the  Southern  Bering  Sea,  Reynolds  et 
al.  (1985)  estimated  that  ice  floes  drift  to 
the  right  of  the  wind  by  approximately  30°  at 
about  4%  of  the  wind  speed  at  3m.  The  ice 
breeze  is  analogous  to  the  landbreeze.  There¬ 
fore,  1°K/ 100km  surface  temperature  gradient 
roughly  produces  0.5  m/s  ice  breeze,  which  in 
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turn  generates  about  2  cm/s  ice  flow.  Pa¬ 
quette  and  Bourke  (1979),  during  a  summer 
cruise  in  the  Chukchi  Sea,  found  the  sea  sur¬ 
face  temperature  to  change  by  over  6°K  across 
a  distance  of  about  25  km  as  the  ice  concen¬ 
tration  fell  from  seven-eighths  cover  to  zero. 
This  is  an  enormous  gradient  by  most  oceanic 
standards  (McPhee,  1983).  Under  such  a  condi¬ 
tion  it  is  found  that  the  ice  breeze  is  around 
12  m/s  (Schmidt  's  estimation),  and  the  asso¬ 
ciated  ice  drift  is  about  48  cm/s  (Reynold  et 
al  's  estimation).  The  coupled  air-ice-ocean 
model  contains  three  parts:  a  thermally 
forced  boundary  layer  air  flow,  a  mechanically 
driven  ice  drift,  and  a  reduced  gravity  ocean. 
The  three  components  are  linked  through  the 
surface  temperature  gradient  and  various  in¬ 
terfacial  stresses. 


THERMALLY  FORCED  BOUNDARY  LAYER  AIR  FLOW 


A  K-theory  approach  planetary  boundary 
layer  air  model  with  modified  Boussinesq  ap¬ 
proximation  and  with  constant  stress  sublayer 
treated  by  Kuo  (  1973)  and  Chu  (  1986, 
1987a, b,c)  is  employed  to  compute  surface  wind 
(ice  breeze'  thex-mally  driven  by  surface  temp¬ 
erature  gradient  across  the  ice  edge.  The  x- 
axis  is  in  the  cross  ice-edge  direction  (posi¬ 
tive  iceward),  and  the  y-axis  is  parallel  to 
the  ice  edge,  as  shown  in  Fig.  1.  It  is  con¬ 
sidered  that  spatial  variations  in  the  MIZ  are 
much  large  perpendicular  to  the  ice  edge  than 
parallel  to  it,  and  hence  derivatives  with  re¬ 
spect  to  y  are  assumed  to  be  zero.  As  dis¬ 
cussed  by  Chu  ( 1986)  the  waterward/iceward  mi¬ 
gration  of  the  MIZ  increases/decreases  the 
surface  temperature  gradient.  The  effects  of 
ice  flow  on  the  surface  temperature  gradient 
can  be  depicted  by  (Chu,  1987b) 

fa'/dtdx  =  (itDTQ/L2)u(x,t)  (1) 

where  0^  is  a  surface  temperature,  u(  is  an  ice 
drift  velocity  in  the  cross  ice-edge  direc¬ 
tion,  DT„  is  a  characteristic  surface  tempera¬ 
ture  differnce  across  the  MIZ,  and  L  is  a 
length-scale  twice  the  MIZ  width.  Subscripts 
'a'  and  * i '  denote  atmosphere  and  ice,  respec¬ 
tively.  The  surface  temperature  is  separated 
into  a  steady  state  and  a  time-dependent  part. 
Both  parts  are  decomposed  into  Fourier  sine 
series.  The  fact  that  the  waterward  monotoni- 
cally  increasing  surface  temperature  often  ob¬ 
served  near  the  MIZ  (Muench,  1983)  implies 
that  the  steady  state  only  has  the  first  mode, 
i.  e.  , 


0a(  x,  t)=^0ak(  t)sin  k*x/L  ♦0asin  nx/b  (2) 

IT 

where  0jk  and  0a  are  Fourier  coefficients.  In¬ 
tegrating  Kuo's  (1973)  planetary  boundary  lay¬ 
er  model  with  slip  and  kinematic  boundary  con¬ 
ditions  at  the  surface  and  finite  condition  at 
the  top,  the  thermally  driven  surface  ’wind  is 
given  by  (Chu,  1986,  1987a, b,c) 


u.(x,t)=VUak;ik(t) 


cos  knx/L  +ualcos  nx/b 


(3) 


va(  x,t)=^vBk;jk(  t)cos  knx/L  +valcos  nx/b  (4) 

where  (u.i(v0)  is  a  surface  wind  vector,  and  uik 


and  vak  are  constants  defined  by 

u*  s  U.J  Vt,; 

vai  E  u< bk-2fu  V  akj/ilii).  k=l , 2, . . .  (5) 

3=1 

where  fu=sin<|>,  tp  is  the  latitude,  and 

U  *  goDTpAZUM^)  (6) 

is  the  scale  of  ice  breeze  wind.  Here 
8=(  v/!l)1//2  is  the  Ekman  depth,  v  is  the  air 
vertical  eddy  viscosity,  and  0i(l  is  a  charac¬ 
teristic  air  temperature  at  surface.  Accord¬ 
ing  to  Kuo  (1973),  v  is  taken  as  5  nr/s.  The 
eigenvalues  /.kj  are  the  roots  with  negative 
real  parts  of  following  sixth  order  algebraic 
equations: 

>.,,+4f(|2>.2  -4k2it2Ri=0,  k=l ,  2 _  (7) 

where 


Ri  e  (  6N/2LSI) 2 


(6) 


is  the  air  Richardson  number.  N  is  the  air 
Brunt-Vaisaila  frequency.  The  integral  con¬ 
stants  ak-  and  bk  in  (5)  are  the  roots  of  the 
following  nonhomogeneous  linear  algebraic 
equations 


=  0, 


i  v 1  -  **«> 


fii 


=  0. 


-2fo  I  ivV<i-*w>  *  bL  =  °< 

3=1 


2knRi  V  akj/>1tj  =  -1/k. 

3-1 


(9) 


where  M  is  a  measure  of  the  effective  depth  of 
the  constant  stress-sublayer  (Kuo,  1973): 
M=v/(CaiU5),  here  Cnj  is  an  air  drag  coeffient 
on  ice.  The  ice-ocean  system  is  mechanically 
driven  by  the  surface  wind,  therefore,  ice  ve¬ 
locity  should  have  the  same  Fourier  components 
as  the  wind,  i.  e.  , 


Uj(  x,  t)=^ujk(  t)cos  knx/L+ucos  nx/b  (10) 


Vj(  x ,  t )  =^]vik(  t )  cos  kux/L+vcos  nx/b  (11) 

where  (ir.v.)  is  the  ice  velocity  vector,  U;., 
vik,  u,  and  v.  are  Fourier  coefficients.  ^(t) 
in  (3)-(4)  is  defined  as  nondiraensionai '  ice 
edge  displacement  (Chu,  1986,  1987b),  i.e. , 


Ld^/dt  =  u^t)  (12) 

If  f()  (i.e..  Latitude),  Ri,  and  M  are  known, 
the  eigenvalues  }.k.  and  the  integral  constants 
a.  and  bk  are  easily  obtained  by  solving  the 
algebraic  equations  (7)  and  (9). 
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rig. 1.  The  coordinate  ayetem  and  air  deflec¬ 
tion  angle  u.  Ice  turning  angle  p 


FREE  ICE  DRIFT  MODEL 

Roed  and  O'Brien  (  1983)  showed  that  the 
internal  ice  stress  and  the  nonlinear  terms 
aren't  crucial  for  the  MIZ  dynamics,  there¬ 
fore,  the  ice  is  nearly  free  drift.  The  mo¬ 
mentum  equations  of  the  ice  drift  are  given  by 

*»/dt  -  fv.  =  TjjU,  +  7»t(  UW-1A )  (13) 

dv./dt  +  fu,  =  7, jV,  ♦  7wl(Vw-V,)  (14) 

where  the  terms  on  the  right-hand  side  repre¬ 
sent  the  water  and  air  stresses  on  ice,  re¬ 
spectively.  Here 

7„  3  l’aCaiU/(  l\tt, )  *  7»i  3  >,«,)  (15) 

Cs  is  a  dimensional  (m/s)  water  drag  coeffi¬ 
cient  on  ice  ( Reod  and  O'Brien,  1983),  and  pn, 
p. ,  and  pu  are  the  densities  of  air,  ice,  and 
water,  respectively.  H  is  a  mean  ice  thick¬ 
ness. 

REDUCED  GRAVITY  OCEAN  MODEL 

Suppose  that  the  ocean  is  composed  of  two 
layers,*  in  which  the  lower  layer  is  deep 
enough  for  motion  in  the  lower  layer  to  be 
vanishingly  small.  Such  a  model  is  referred 
to  a3  a  reduced  gravity  model.  Furthermore, 
the  model  has  the  simplest  form  of  dissipa¬ 
tion,  namely,  Rayleigh  friction  with  decay 
rate  R  and  Newtonian  cooling,  also  with  decay 
rate  R.  The  momentum  and  continuity  equations 
for  the  ocean  are  written  by  (Gill,  1982;  Roed 
and  O'Brien,  1983) 


( 1. e.  ,  the  fraction  of  area  covered  by  ice). 
Hibler  and  Ackley  (1983)  took  the  50%  limit  of 
A  as  the  ice  edge.  Inside  the  MIZ,  the  ice 
concentration,  A,  decreases  waterward  very 
slowly  from  the  boundary  between  pack  ice  and 
the  MIZ  (where  A~l)  to  some  place  near  the 
ice  edge,  and  then  diminishes  very  rapidly  to 
the  ice  edge.  If  the  area  we  focus  on  is  not 
very  close  to  the  ice  edge,  we  may  set  A~l. 
The  portion  of  air  stress  on  the  water, 
(l-A)ywVa,  is  then  neglected  compared  to  ice 
stress  oh  the  water.  Ay.  ( V,-Vw ) .  Setting  A=1 
and  eliminating  hw  from  (16)  and  (16),  the  mo¬ 
mentum  equations  ('  16 )  and  (17)  become 

( d/dt+  R ) uw- f vw- (  g*Hw/RL* )d2uw/dx*=ylw( u.-uw )  ( 19) 

(d/£>t+R)vw+fuw  =7iw(vj-vw)  (20) 

SOLUTIONS 

If  the  ice-ocean  is  considered  as  one 
system,  the  air  stress  is  the  only  forcing 
term  in  (13),  (14),  (19),  and  (20).  The  ice 

velocity  and  the  ocean  current  should  have  the 
same  Fourier  components  as  the  surface  wind 
since  our  system  is  linear.  According  to  (3) 
and  (4)  the  solutions  have  the  following 
forms: 


U.(  x,  t)=Vu.kexp(  pkt)coskitx/L  +ujcosnx/L  (21a) 

X 

v,(  x,  t)=^_vjkexp(  pkt)cosknx/L  +v|cosrx/L  (21b) 
uw(  x,  t)=^]uakexp(  pkt)coskax/L  +uwcosnx/L  (21c) 

vw(  x,  t)=^vwkexp(  pkt)cosknx/L  +vivcositx/L  (21d) 

k  _ 

where  a.,  v|k,  u  k,  v  k  (k=l,2,...),  u,,  v.,  u  , 
and  vw  are  the  Fourier  coeffients,  and  pk 
(k=l,2,...)  are  the  eigenvalues  of  the  sys¬ 
tem.  Substituting  (21a)-(21d)  into  (13), 
(14),  (19),  and  (20)  we  have  a  set  of  nonhomo- 
geneous  linear  algebraic  equations  for  the 
steady-state  Fourier  coefficients  (0w,  vw,  uj( 


and  Vj): 

[R+rt2g*/(RL:)+C|w]Gw  -f?w-7iwGj  =  0,  (22a) 

fQw*(R+Ciw)vw  -7it,v=  0,  (22b) 

-I»iSw+V,G,  -f9i  “  TAi  (22c) 

“ViVvA  +fQ,  =  7aiV#1  (22d) 


and  the  dispersion  relations  for  the  eigenva¬ 
lues  p.: 


( ^t+RJu^.-fv^  =-q*c\./dx+{  1-A)7awua 


*AYiw(u,-u  ) 


( d/dt+R) vw+  fuw  =(  1-A ) 7awva  +A/jw(  vrvw )  (  17 ) 

( <3/i?t+R)hw  +Hvvcuw/5x  =  0  (18) 

where  the  second  and  third  terms  on  the  right- 
hand  side  of  (16)  and  (17)  represent  the  air 
and  ice  stresses  on  water,  and  g*  is  a  reduced 
gravitational  acceleration.  H,v  the  equilibrium 
depth  of  the  pycnocline,  hw  the  thickness  of 
the  upper  layer,  7i,v3Cw,/H,1 ,  and  7awsl'.C;mU/ 
pH,  Cu  i3  an  air  drag  coefficient  on  water. 
a"  is  the  ice  compactness  or  ice  concentration 


-f  ^u^/L  /wi  o 

f  Mk*7wi  7aiVlk/L  0  Yw1 

1  0  |.k  0  0  =0, 

7iw  0  0  Mk+Rj  -f 

0  'hr  0  f  p^ 

where  R,=R+Cm,,  R,=Rt,»g*Hwirk;:/RLJ. 


(23) 
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MEAN  ICE  DIVERGENCE/CONVERGENCE  CRITERION 

In  this  section  we  discuss  the  steady 
state  solutions  (22a)-(22d).  Chu  (1937c)  shows 
that  the  absolute  value  of  ice  breeze  deflec¬ 
tion  angle  (  angle  between  surface  temperature 
gradient  and  surface  wind),  |u|,  increases 
with  increasing  latitude  and  decreases  with 
increasing  Ri.  At  65°  (N  and  S)  latitude  it 
varies  from  41.8"  when  Ri  =  0.001  to  u  =  27.5" 
when  Ri  =  0. 04.  The  angle  between  surface 
wind  and  ice  velocity,  |t  (as  shown  in  Fig.  2), 
is  defined  by 

(I  =  tan"1(v/ui)  -  u  (24) 

where  uf  and  v,  are  the  roots  of  linear  alge¬ 
braic  equations  (22a)-(22d).  An  ice  diver¬ 
gence/convergence  criterion  is  simply 

>  n/2  ice  convergence 

l“+PI  (25) 

<  n/2  ice  divergence 

Fig.  2  indicates  the  mean  ice  divergence/con¬ 
vergence  criterion, 

tan-1(vj/ui)  =  n/2  (26) 

for  different  Ri  in  the  (H,Hw)  plane.  The 
curves  separate  the  parameter"  plane  ( H , Ha) 
into  ice  convergence  and  ice  divergence  pa!rts. 
We  can  differentiate  between  ice  divergence 
and  ice  convergence  due  to  the  ice  breeze  in 
the  MIZ  by  the  parameters  Ri ,  H  ,  and  H  Ice 
divergence  appears  in  the  smal'l  H.  and'  large 
Hw  area,  however  on  the  contrary,  ice  conver¬ 
gence  appears  in  the  large  H-  and  small  H  re¬ 
gion.  The  critical  curve  (  |o+p|=n/2)  moves 
from  lower  left  to  upper  right  in  the  (H(,H.) 
plane  from  less  to  more  stable  atmosphere.'  If 
the  drag  coefficient  C  .  is  doubled,  the  turn¬ 
ing  angle  p  only  has  minor  changes  ( less  than 
10%). 


R  f  Ml 

V 

Fig.  2.  H-Hw  curvet  for  zero  ice  divergence/ 
convergence  for  five  different  values  of  Ri 
< after  Chu.  1987a). 


TWO  BIFURCATIONS  OF  TIME-DEPENDENT  ICE  DRIFT 

Ice  drift  observations  in  the  Greenland 
Sea  from  1978  April  28  to  September  3  (Vinje, 
1982)  show  two  different  types  of  ice  motion: 
oscillatory  and  nonoscillatory.  The  mesoscale 
air-ice-ocean  feedback  mechanism  may  explain 
this  phenomena.  If  the  ocean  is  deep  enough 
such  that  Clw  can  be  neglected  against  R,  the 
dispersion  relation  (23)  is  simplified  as 

^+R+itVg*Hw/RL2  -f 

=  0,  (27) 

f  Mk+R 

representing  inertial-gravity  water  waves,  and 

“k+Cwi  -f  CaiUk«/L 

£  *\+Cwi  Caiv,a/L  =0  (28) 

1  0  Mk 

showing  the  ice  drift  patterns.  The  roots  nk 
(k=l,2,...)  of  the  cubic  equations  (28)  give 
the  e-folding  time  dependence  of  the  kth  com¬ 
ponent  of  ice  velocity.  The  instability  cri¬ 
terion  for  the  kth  mode  of  ice  motion  in  the 
MIZ  can  be  written 

Re(pk)  <  0  (29a) 

the  kth  mode  of  ice  velocity  decreases  with 
time, 

Re(,ik)  =  0  (29b) 

neutral,  and 

Re(nk)  >0  (29c) 

the  kth  mode  of  ice  velocity  increases  with 
time.  We  define  the  time-increasing  mode  of 
ice  velocity  as  an  unstable  mode. 

The  oscillation  criterion  for  the  kth 
mode  of  ice  motion  is  given  by 

In*( **k >  =  (nonoscillatory)  (30a) 

Im(fik)  *  0,  (oscillatory)  (30b) 

We  compute  all  the  roots  of  (28)  for  different 
values  of  the  parameters  K,  N,  DTn,  and  obtain 
three  roots  M),  H?,  |*3  at  each  point  of  the  pa¬ 
rameter  space  (R(,N-",  DTn).  Here  M  is  real 
throughout  that  space,  and  n2  and  ji,  are  real 
for  some  values  of  (H  N:,DT„)  and  are  complex 
else  where. 

Fig.  3  shows  the  surface  |i.=0  for  k=l  in 
the  three  dimensional  space  (R,N:,DT(|).  This 
surface  divides  the  space  into  two  parts  cor¬ 
responding  to  growing  and  decaying  modes.  The 
growing  mode  generally  appears  when  DT(|  ex¬ 
ceeds  some  critical  value  around  10°C.  This 
critical  value  decreases  with  R,  and  becomes 
very  small  (around  1°C)  in  the  region  of  small 
mean  ice-thickness  (H.<1.5  m).  Ice  motion 
corresponding  to  the  eigenvalue  p  is  nonos¬ 
cillatory. 
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FX9.  3.  Separation  of  nonoseiiiatory  decaying 
and  growing  modes  referring  to  real  eigenva¬ 
lues  a,  for  *=1  (after  Chu.  1986). 


Ilf,  ®'F!r,:lon  o£  nonoseiiiatory  and 


and 


o«- 

i*: 


Fig.  4  indicates  the  surface  of  Re(p,)=0 
(or  Re(|i,)=0)  for  k=l  in  the  parameter  space 
(H,  N2,DT,j).  This  surface  separates  the  whole 
space  into  growing  and  decaying  parts.  The 
growing  mode  is  present  when  DTn  exceeds  some 
critical  value  which  .is  a  function  of  H,  and  N* 
(when  E=  4.  5  m  and  N:=0. 00011/s*,  the  critical 
value  is  3°C). 


Fig. 4.  Separation  of  decaying  and  growing 
nodes  referring  to  eigenvalues  ji.  and  ny 


Fig.  5  reveals  the  segregation  of  nonos- 
cillatory  and  oscillatory  modes  relating  to  ji, 
and  jiv  Comparing  Fig.  5  with  4,  we  find  that 
the  decaying  mode  of  p.,  and  p,  is  generally  no- 
noscillatory  whereas  ‘the  growing  mode  of  p, 
and  p,  is  generally  oscillatory. 


Fig. 3-5  show  the  following  results:  (i) 
both  nonoseiiiatory  and  oscillatory  modes 
share  a  common  area  ( restricted  to  the  region 
of  small  DT„)  in  the  parameter  space;  however, 
(ii)  nonoseiiiatory  and  oscillatory  growing 
modes  occupy  different  regions  in  the  parame¬ 
ter  space.  The  segregation  of  the  two  modes 
depends  largely  on  DT()  and  H. 

Whether  ice  motion  grows  or  decays  sub¬ 
stantially  depends  on  a  first  critical  yalue 
of  the  parameter  DTn  (when  N=l. 45xl0”2s”i  and 
a=2.  5m,  this  critical  value  is  S‘’C).  When  DT„ 
is  smaller  than  the  first  critical  value,  the 
motive  force  is  so  small  that  it  cannot  over¬ 
come  dissipative  effects  and  does  not  make  ice 
motion  unstable.  If  DT„  becomes  large  enough 
to  overcome  the  dissipative  effects  of  fric¬ 
tion,  ice  motion  becomes  larger. 

Whether  ice  motion  is  oscillatory  or  no- 
noscillatory  largely  depends  on  DT(I  and  the 
properties  of  ice.  If  DTU  exceeds  the  first 
critical  value  but  does  not  reach  a  second 
critical  value  which  mostly  depends  on  N 
(i.e. ,  when  M=0. 0145/s,  the  second  critical 
value  is  14“C),  and  when  ice  is  thin  (general¬ 
ly  during  summer)  the  ice  motion  is  nonoscil- 
latory,  however,  when  the  ice  is  thick  (gener¬ 
ally  during  winter)  the  ice  motion  is 
oscillatory.  If  DT(.  exceeds  the  second  criti¬ 
cal  value,  only  the  nonoseiiiatory  growing 
mode  appears. 

For  the  nonoseiiiatory  growing  mode  the 
time  during  which  the  ice  doubles  its  speed  is 

T,  =  ln(  2/p, ) ,  (  p,  >0 )  (31) 

The  doubling  time  treated  as  a  function  of  DTI( 
(for  Fl=lm,  M=0. 01/s)  is  shown  in  Fig.  6,  which 
displays  a  decrease  of  T|  with  an  increase  of 
DTn.  T,  changes  from  2.  2  to  0.  18  hr  as  DTU  var¬ 
ies  from  5°  to  22"C. 
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ICE  EDGE  UPWELLING 

Buckley  et  al.  (1979)  reported  a  pro¬ 
nounced  wind-driven  upwelling  along  the  edge 
of  the  ice  pack  in  an  expedition  into  the  Arc¬ 
tic  Ocean  north  of  Spitsbergen.  The  aim  of 
this  section  is  to  investigate  the  effect  of 
the  air-ice-ocean  feedback  mechanism  on  ice 
edge  upwelling.  If  the  sinusoidal  type  surface 
temperature  perturbation  across  the  MIZ 
doesn  t  change  as  the  ice  edge  moves  in  a 
speed  uj(  the  parameterization  of  surface 
thermal  condition  (1)  becomes 

»n(x.t)=OaU  -DT0sin*[  (  x-Ujt)/L]  (33) 

The  negative  sign  in  (34)  comes  from  the  x- 
axis  pointing  iceward  (Fig.  1).  Let 
c=(x-Ujt)/L(  the  surface  wind  driven  by  the 
horizontal  temperature  gradient  (33)  is  calcu¬ 
lated  by 


Fig.  6  Dependence  of  doubling  tin*  on  DT.  for 
the  npnoocillat.ory  growing  mod*  for  H=lm  and 
N=!0  *  (after  Chu,  1986).  ' 

The  oscillatory  growing  mode  is  ice  mo¬ 
tion  due  to  m  (or  |i,),  whose  real  part  is  po¬ 
sitive  and  the  imaginary  part  is  different 
from  zero.  The  growing  oscillatory  mode  is  lo¬ 
cated  in  the  area  where  DT0  is  larger  than  the 
first  critical  value  and  smaller  than  the  sec¬ 
ond  critical  value.  As  in  the  nonoscillatory 
case  the  growth  rate  Re(|i.,/})  increase  with  an 

increase  of  DT0  (increase  of  forcing).  The 
period  of  the  oscillatory  mode  is  defined  by 

Tp  =  2n/|Im(M:,3)|  (32) 

Fig.  7  shows  the  doubling  time  and  the  period 
as  a  function  of  DT0  for  fi=5m  and  N=0. 01/s. 
The  doubling  time  decreases  monotonical ly  with 
an  increase  of  DT„.  However,  the  period  in¬ 
creases  slightly  when  DT(|  varies  from  5°C  to 
10°C,  and  then  increases  rapidly  with  DT(). 
When  DT0  =  12.5°C  the  period  is  nearly  one 
day. 


Dy6c) 


ua!z=0  =  Ualcos  *4.  (34) 

Valz=0  =  Va|cos  "4,  (35) 

where  ua|  and  v  .  are  computed  by  (5).  Neg¬ 
lecting  inertial  oscillation  and  assuming  the 
ocean  velocities  to  be  small  ( Roed  and  O'¬ 
Brien,  1983),  the  momentum  equations  for  the 
reduced  gravity  ocean  model  (19)  and  (20)  are 
written  by  '  1 


Ruw-*vw=-g*ava*+l  l-H<  4)  1  C„ualcos«4*Ciw(  u,-u  )H(  4 ) 

(36) 


( <5/f’t+R)vw+fuw-[  1-H(  4)lCawva|cos*4+Cjw(  v-v  )H(£) 

(  37) 

where  H(£)  is  the  Heaviside  function  defined 

as 


1,  if  4>0,  (MIZ) 

H(4)  (38) 

,  0,  if  £<0,  (open  water) 

The  equations  for  the  free  drift  ice  (13), 
(14)  and  for  the  reduced  gravity  ocean  (36), 
(37),  and  (18)  are  the  basic  equations  for 
this  section.  The  solution  is 


-(Hw//.)[C1exp(-LI4IA) 

+  (>.n/L)(  p(|cos)t4— Pusina4)]  ,  (4<0) 

dhw/Bt  (39) 

,  - (  Hw/>. )  (  C,exp(  - L  1 4  |  />. ) 

+  (>.it/L)(p,1cosn4-p,,sinnd))  ,  (4>0) 

where  the  parameters  are  defined  by 

*  =  )1/2/f, 

P)  |  =  *CWU,U„  /  f  2i.t  l*  (  1*/L  )  2 1  . 

P|2  =  cawval/f[l+(>.n/L)2)  , 

P2)  =  nC.wU,J/f2M  l  +  (  >.n/L  )2]  , 

Pj2  =  ciwV/ f  t  l+(  ln/L)2]  , 

ci  =  <P22-Pij)/2  +  >"(P:i-Pm)/(2L) 


Fig. 7.  Dependence  of  doubling  time  and  period 
on  DT,  fox  qsciiiatory  growing  mode  for  K=5m 
and  Ns10"fca"*  (after  Chu.  1986). 
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(40) 


c,  =  -<P22-Pl2)/2  +  '•*<  P’i  ~P)  i )/( 2L) 

where  Uj,  Vj  are  computed  from  nonhomogeneous 
linear  algebraic  equations  (22a)-(22d). 
Eq. (39)  shows  that  the  ice  breeze  disturbance, 
which  is  local  to  the  MIZ,  will  cause  the  re¬ 
duced  gravity  ocean  model  to  adjust  to  differ¬ 
ential  Ekman  transport  (both  from  changes  in 
stress  at  the  ice/ocean  boundary  and  from  hor¬ 
izontal  variation  in  wind)  by  increasing  or 
decreasing  thickness.  The  rate  of  shallowing 
(thickening)  is  equivalent  to  an  upwelling 
( downwelling) .  This  effect  of  ice  breeze  on 
ice-edge  upwelling/downwelling  can  be  parame¬ 
terized  in  terms  of  a  differential  surface 
heating  condition,  i.e.  , 

^hw/c’t=A,exp(  -LI4I  A)+a2[  £ea/^-B(  0a-eaQ)]  (  41) 
where 

A,  =  -(HWC,A)I  l-H(i)]  ♦  (  HWCJA)H( ;)  , 

A,  =  (Hw/LDT0)r,,l  1-HU)]  ♦  fhjHU),  (42) 

B  =  *(P|3/P|i>[l~H(t)]  *  (P22/P2I>H<^- 

Eq.  (41)  shows  the  relationship  between  ice- 
edge  upwelling/downwelling  and  the  surface 
temperature  distribution  across  the  ice  edge. 

We  compute  dh^y^t  as  a  function  of  x  and 
t  for  H.=l,  and  6m,  respectively,  as  shown  in 
Fig.  8.  'it  is  seen  that  upwelling  ( <?hw/<Jt  <  0) 
appears  near  the  ice  edge.  The  minimum  value 
of  dh^/dt  (maximum  upwelling)  reaches  18.3  m 
day-1  for  H.  =  lm,  and  19.8  m  day-1  for  H.=6  m. 
If  the  distance  between  two  lines  of  -5  m 
day"1  is  taken  as  the  width  of  ice-edge  up¬ 
welling,  it  is  found  that  this  width  decreases 
with  the  increase  of  mean  ice-thickness 
The  upwelling  width  is  around  40  km  for  Eg  = 
lm,  and  25  km  for  H  =  6m.  These  values  are 
quite  comparable  to  the  baroclinic  radius  of 
deformation  i.  ( 23  km). 


ACROSS  EDGE  DISTANCE  (KM  ) 


Fig.  8.  hi /it  (m  as  a  function  of  x 

and  t  for  (a)  H.,  =  l  m,  and  ( b )  H  =  6  m. 


CONCLUSIONS 

(  a)  This  air-ice-ocean  coupled  model  is 
intended  to  depict  only  the  mesoscale  process¬ 
es  of  air-ice-ocean  interactions  in  the  MIZ. 
The  synoptic  scale  pressure  gradient  may  addi¬ 
tionally  produce  surface  winds  in  the  MIZ,  and 
large-scale  ocean  current  may  also  drive  ice 
drift.  These  processes  are,  however,  beyond 
the  scope  of  this  paper.  Nevertheless,  where 
the  ice-to-open-water  temperature  gradient  is 
strong,  the  mesoscale  feedback  mechanism  dis¬ 
cussed  here  may  become  as  strong  as,  or 
stronger  than,  the  synoptic  scale  and  oceanic 
forcings. 

(b)  The  model  shows  the  different  effects 
of  ice  breeze  on  ice  flow,  and  gives  out  an 
ice  divergence/convergence  criterion  in  the 
MIZ. 

(c)  The  ice  motion  has  two  bifurcations. 
First,  it  bifurcates  into  a  decaying  or  grow¬ 
ing  mode,  which  depends  in  most  cases  on  the 
mean  surface  temperature  difference  DT„  repre¬ 
senting  the  strength  of  the  forcing.  When  DT() 
is  small,  the  decaying  mode  exists.  However, 
when  DT0  exceeds  the  first  critical  value,  the 
growing  mode  appears.  Secondly,  the  growing 
mode  bifurcates  into  nonoscillatory  and  oscil¬ 
latory  states  depending  on  DT„  and  the  proper¬ 
ties  of  ice. 

(d)  Surface  winds,  thermally  generated  by 
the  temperature  gradient  across  the  ice  edge, 
cause  the  reduced  gravity  ocean  model  to  ad¬ 
just  to  differential  Ekman  transport  both  from 
changes  in  stress  at  the  ice/ocean  boundary 
and  from  horizontal  variation  in  winds.  Such 
differential  Ekman  transport  generates  the  ice 
edge  upwelling. 
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